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Abstract. 

In the framework of the theory of open systems, we give a description of quantum 
entanglement and quantum discord for two non-interacting modes embedded in a 
thermal environment. We describe the evolution of entanglement in terms of the 
covariance matrix for Gaussian input states. For all values of the temperature of the 
thermal reservoir, an initial separable squeezed thermal state remains separable for 
all times. In the case of an entangled initial squeezed thermal state, entanglement 
suppression (entanglement sudden death) takes place, for non-zero temperatures of 
the environment. The Gaussian quantum discord, which is a measure of all quantum 
correlations in the state, including entanglement, decays asymptotically in time. 

PACS numbers: 03.65.Yz, 03.67.Bg, 03.67.Mn 

1. Introduction 

In the framework of the theory of open systems based on completely positive quantum 
dynamical semigroups, we study the dynamics of quantum entanglement and quantum 
discord of a subsystem consisting of two uncoupled modes interacting with a common 
thermal bath. The initial state of the subsystem is taken of Gaussian form and 
the evolution under the dynamical semigroup assures the preservation in time of the 
Gaussian form of the state. In Sec. 2 we give the solution of the evolution equation 
for the covariance matrix corresponding to a generic two-mode Gaussian state of the 
two uncoupled modes interacting with the environment. We investigate in Sec. 3 
the dynamics of entanglement, by analyzing the time evolution of the logarithmic 
negativity, which characterizes the degree of entanglement of the state. We study also 
the time evolution of the Gaussian quantum discord, which is a measure of all quantum 
correlations in the bipartite state including entanglement. Summary is given in Sec. 4. 

2. Equations of motion for two modes interacting with an environment 

We study the dynamics of a subsystem composed of two non-interacting modes in weak 
interaction with a thermal environment. In the axiomatic formalism based on completely 
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positive quantum dynamical semigroups, the Markovian irreversible time evolution of 
an open system is described by the Lindblad-Kossakowski master equation. We are 
interested in the set of Gaussian states, therefore we introduce such quantum dynamical 
semigroups that preserve this set during time evolution of the system. 

The two-mode Gaussian state is entirely specified by its covariance matrix, which 
is a real, symmetric and positive matrix: 
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where A, B and C are 2x2 Hermitian matrices. A and B denote the symmetric 
covariance matrices for the individual reduced one-mode states, while the matrix C 
contains the cross-correlations between modes. The equations of motion for the quantum 
correlations of the coordinates x, y and momenta p x , p y of the two modes are the following 
(T denotes the transposed matrix) pQ: 
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m is the mass, u% and U2 are the frequencies of the non-resonant modes, and diffusion 
coefficients D xx , D xpx ,... and dissipation constant A are real quantities. 
The time-dependent solution of Eq. (|2]) is given by [1] 

a(t) = M(t)[a(0) - o-(oo)]M T (t) + a(oo), (5) 

where the matrix M(t) = exp(Yt) has to fulfill the condition lim^oo M(t) = 0. The 
values at infinity are obtained from the equation 

Ya(oo) + a(oo)Y T = -2D. (6) 



3. Dynamics of continuous variable entanglement and discord 

3.1. Time evolution of entanglement and logarithmic negativity 

In order to quantify the degree of entanglement of the two-mode states it is suitable 
to use the logarithmic negativity. For a Gaussian density operator, the logarithmic 
negativity is given by = max{0, — log 2 2z>_}, where z>_ is the smallest of the two 



Entanglement and discord in open quantum systems 



3 



symplectic eigenvalues of the partial transpose a of the 2-mode covariance matrix a. In 
our model, the logarithmic negativity is calculated as [21 El IH [5] 
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We assume that the initial Gaussian state is the two-mode squeezed thermal state, 
with the covariance matrix of the form [6] 
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with the matrix elements given by 



a = rii cosh r + U2 sinh r + - cosh 2r, 

b = n\ sinh 2 r + n 2 cosh 2 r + - cosh 2r, 

1 . . 
c = -yrii + ^2 + 1) sinh2r, 



(10) 
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where n\, n 2 are the average numbers of thermal photons associated with the two modes 
and r denotes the squeezing parameter. In the particular case n\ — and n 2 = 0, 
(|9| becomes the covariance matrix of the two-mode squeezed vacuum state. A two- 
mode squeezed thermal state is entangled when the squeezing parameter r satisfies the 
inequality r > r s [6], where 

(m + l)(n 2 + l) 



cosh r s 



(13) 



ni + ri2 + 1 

We suppose that the asymptotic state of the considered open system is a Gibbs 
state corresponding to two independent modes in thermal equilibrium at temperature 
T. Then the quantum diffusion coefficients have the following form (we put h = 1) [7]: 
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The calculations with separable initial squeezed thermal states show that these 
states remain separable for all values of the temperature T and for all times. The 
evolution of entangled initial squeezed thermal states with the covariance matrix given 
by Eq. ^ is illustrated in Figure 1, where we represent the dependence of the 
logarithmic negativity E^(t) on time t and temperature T for the case of an initial 
symmetric Gaussian state, when a = b. We observe that for a non-zero temperature 
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Figure 1. Logarithmic negativity versus time t and temperature T for an 
entangled initial symmetric squeezed thermal state with squeezing parameter r = 2, 
rii = l,ri2 = 1 and A = 0.1, Wi = W2 = 1- We take m = h = k = 1. 



T, at certain finite moment of time, which depends on T, E^it) becomes zero and 
therefore the state becomes separable. This is the so-called phenomenon of entanglement 
sudden death. It is in contrast to the quantum decoherence, during which the loss of 
quantum coherence is usually gradual [HI E]- For T = 0, E N (t) remains strictly positive 
for finite times and tends asymptotically to for t — > oo. Therefore, only for zero 
temperature of the thermal bath the initial entangled state remains entangled for finite 
times and this state tends asymptotically to a separable one for infinitely large time. 
Dissipation favorizes the phenomenon of entanglement sudden death - with increasing 
the dissipation parameter A, the entanglement suppression happens earlier. The same 
qualitative behaviour of the time evolution of entanglement was obtained previously 
[TUt [TT] in the particular case n\ — and n 2 = corresponding to an initial two-mode 
squeezed vacuum state. 

3.2. Gaussian quantum discord 

Quantum discord was introduced [12] as a measure of all quantum correlations in 
a bipartite state, including - but not restricted to - entanglement. Originally the 
quantum discord was defined and evaluated mainly for finite dimensional systems. Very 
recently p31 [H] the notion of discord has been extended to the domain of continuous 
variable systems, in particular to the analysis of bipartite systems described by two- 
mode Gaussian states. Closed formulas have been derived for bipartite thermal squeezed 
states [13] and for all two-mode Gaussian states [H]. The Gaussian quantum discord 
of a general two-mode Gaussian state can be defined as the quantum discord where 
the conditional entropy is restricted to generalized Gaussian positive operator valued 
measurements (POVM) on the mode 2, and in terms of symplectic invariants it is given 
by (the symmetry between the two modes 1 and 2 is broken) [2] 
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Figure 2. (a) Gaussian quantum discord D versus time t and temperature T for an 
entangled initial symmetric squeezed thermal state with squeezing parameter r = 2, 
rii = ri2 = 1 and A = 0.1, Wi = a>2 = 1. We take m = % = k = 1. (b) Same as in (a), 
with the explicit thresholds shown for 13 = and D = 1. 
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and v T are the symplectic eigenvalues of the state, given by 
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where A = det A + det 5 + 2 det C. Notice that Gaussian quantum discord only depends 
on | det C\, i.e., entangled (det C < 0) and separable states are treated on equal footing. 

The evolution of the Gaussian quantum discord D is illustrated in Figure 2, where 
we represent the dependence of D on time t and temperature T for an entangled initial 
symmetric Gaussian state, taken of the form of a two-mode squeezed thermal state ([9]), 
for such values of the parameters which satisfy for all times the first condition in formula 



(17). For entangled initial states the Gaussian discord remains strictly positive in time 



and in the limit of infinite time it tends asymptotically to zero, corresponding to the 
thermal product (separable) state, with no correlations at all. The Gaussian discord 
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has nonzero values for all finite times and this fact certifies the existence of nonclassical 
correlations in two-mode Gaussian states, either separable or entangled. 

From Figure 2 we notice that, in concordance with the general properties of the 
Gaussian quantum discord [14J, the states can be either separable or entangled for D < 1 
and all the states above the threshold D = 1 are entangled. We also notice that the 
decay of quantum discord is stronger when the temperature T is increasing. 

4. Summary 

In the framework of the theory of open systems we investigated the Markovian dynamics 
of quantum correlations for a subsystem composed of two non-interacting modes 
embedded in a thermal bath. We have described the time evolution of the logarithmic 
negativity in terms of the covariance matrix for squeezed thermal states, for the case 
when the asymptotic state is a Gibbs state corresponding to two independent quantum 
harmonic oscillators in thermal equilibrium. For all values of the temperature of the 
thermal reservoir, an initial separable Gaussian state remains separable for all times. 
In the case of an entangled initial squeezed thermal state, entanglement sudden death 
takes place for non-zero temperatures of the environment. Only for a zero temperature 
of the thermal bath the initial entangled state remains entangled for finite times, but 
in the limit of infinite time it evolves asymptotically to an equilibrium state which 
is always separable. The Gaussian quantum discord decreases asymptotically in time. 
This is contrast with the sudden death of entanglement. The time evolution of quantum 
discord is very similar to that of the entanglement before the sudden suppression of 
the entanglement. After the sudden death of the entanglement, the nonzero values of 
quantum discord quantify the nonclassical correlations for separable mixed states. 
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